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eponsonhlp  of  the  U.S  Department  of  Arey 
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EOTCa  HLADE  VIBRttIZCN  MODES  ;j©  TODESUaiOIBB 


SUUUARI 


Tlia  primary  object  of  the  research  has  boon  the  numerioal 
evaluation  of  certain  exact  solutions  of  the  differential  equation  of 
notion  of  a  flexible  rotor  blade  under  the  notion  of  centrifugal  force 
and  inertia,  giving  the  natural  frequencies  of  free  vibration  and.  the 
defleotion  curves  or  nodes  of  the  different  possible  typeB  of  vibration. 

These  nunerionl  results  apply  to  various  degrees  and  kindn  of  blade 

taper,  and  they  have  been  tabulated  and  plotted.  The  taper  range  is 

(E0)  Uniform  moss  distribution  along  the  blade 
radius 

(Do)  Linear  (triangular)  taper 

(A)  Elliptical  taper 

(B)  Generalised  elliptical  taper 

(c)  Parch  olio  taper 

'D)  Generalised  parabolic  taper 

(i<)  Blade  extension,  parabolic  type  taper 

(E1 )  hyperbolic  taper 

(So)  can  be  regarded  as  a  special  oase  of  (B)  or  (D)  and  (D0)  of  (D), 
while  ()>)  is  s.lnply  (D)  oxtended  across  the  axis  of  rotation.  (E*)  is 
the  Uniting  oase  of  (E) . 

Pron  tfao  results  for  the  mathematical  taper  forms  listed,  the 
frequencies  of  the  various  modes  for  non- mathematical  taper  farms,  as 
used  in  aoet  actual  blades,  oon  be  interpolated  fairly  accurately 
according  to  the  degree  of  taper,  specified  by  the  distance  of  the 
blade  oentrodd  froa  the  tads,  in  terns  of  the  tip  radius.  In  inter¬ 
polating  nodes  (defleotlone)  none  regard  should  also  be  paid  to  the 
type  of  taper,  e  g.  whether  oonoave  or  ooaraac. 

HUde  root  oondltlona  influenoe  the  nodes  and  fireqpenoles. 

In  general  these  oan  be  oorered  by  the  location  of  the  flapping  hinge, 
or  the  equivalent  hinge  ponitlon  for  built-in  blade  roots.  Hinge 
looatlon  is  one  of  the  variables  included  in  the  eslmlatlons,  so 
tint  its  effect,  which  is  sot  very  sensitive  to  taper,  oen  easily  be 
found. 
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For  actual  blades  pm— going  flexural  rigidity  a  itU-faua 
formula  for  rlhration  frequency  laO*  ■  (k*w*  +  l®1)  vbsrs  0  le  the 
actual  resultant  frequency,  (tar)  the  flexible -blade  frequency  considered 
In  this  report,  and  Iqj  the  frequency  of  the  non-rotating  blade  vibrating 
as  a  bean.  The  present  work  shews  that  the  formula  is  exact  for  a 
certain  M  standard"  distribution  of  bend  rigidity  (Ej)  along  the  blade 
radius ,  the  bending  node  then  being  identical  with  the  corresponding 
sods  for  a  flexible  rotating  blade  having  the  ease  Base  taper.  Even 
when  not  exaot  the  formula  represents  a  useful  analysis  of  the  frequency 
shot  zing  the  results  of  a  change  in  r.p.m.  or  flexural  rigidity,  and 
checking  more  da  tolled  calculations. 

Arising  from  the  present  work,  although  not  port  of  the 
original  programs  end  therefore  not  dealt  with  in  any  detail  in  the 
report: 

(i)  The  response  of  a  rotating  blade  to  periodic,  sudden 

or,  as  a  special  case,  static  loading  oan  be  calculated 
exactly  for  any  of  the  tapered,  blades  having  the  "standard" 
stiffness  distribution.  This  provides  a  method  of 
stressing  the  blade. 

(ii)  hltern&tively,  independent  eaaot  solutions  are 

available  for  the  static  loading  case,  oognate  with  those 
for  the  vibrating  blade. 

(ill)  Applicable  not  only  to  the  taper  forms  specified 

above ,  but  also  to  the  wider  taper  range  obtained  by 
truncating  these  forma,  complete  and  exact  solutions 
are  available  for  vibrating  beams  having  a  wide  range 
of  mass  and  atiffnese  taper. 

(iv)  Shear  deflection  effects  oan  readily  be  included 

in  the  calculations  either  far  the  non-rotating  beams 
or  the  rotating  blades. 

(v)  Truncated  uniform  or  tapered  blades  carrying  tip 

susses  are  oovered  by  a  fairly  slapl*  extension  of  the 
solution  for  the  intaot  blade. 

(vi)  The  differential  equation  for  a  rotor  blade, 

idwthsr  statically  loaded  or  vibrating  is  eathematloally 
Identical  with  that  for  a  so-nailed  box  wing  consisting 
of  too  spars  and  a  torsion  box.  Solutions  and  results 
for  such  wings  oan  therefore  be  deduoed  from  those  given 
hire. 
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LIST  01  SYMBOLS 


x  Blade  station  (fraction  of  tip 

radius,  B) 

y  Blade  defleotion 

w  Angular  velocity  of  rotor 

p  Blade  frequency 

k  Blade  frequency  ratio  p/w 

n,  n  indices 

g  Acceleration  due  to  gravity 

P  Centrifugal  tension  in  blade 

ft  Coefficient 

KI  Blade  bend  rigidity 

GA  Blade  shear  rigidity 


i 


Reference! 


Airoraft  Engineering,  Novenbor  1958 
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StnWSB  Off  EKKJUENC2ES 


Values  of  lc*  .  p*/w* 


Maas  Mode 


Case 

Distribution 

Centroid 

0 

_  is 

1ft 

2ft 

Jft 

Bo 

Uniform 

0-5 

1 

3 

6 

10 

15 

21 

28 

E 

(1  -  *•)* 

0-4.56 

1 

2.8 

5-4- 

8-8 

13 

18 

23-8 

E 

(1  -  *)-* 

0-455 

1 

5  38 

13-3 

A 

(1  -  x*)* 

o  .424. 

1 

2-67 

5 

8 

11-67 

16 

21 

E 

i 

rH 

s_/ 

0-404. 

1 

a. 

5  05 

12-2 

C 

(l  -  *)* 

0-40 

1 

4-- 75 

11-07 

B 

(1  -  X*) 

0-375 

1 

2-5 

4-5 

7 

10 

13-5 

17-5 

Do 

(1  -x) 

0-333 

1 

4-2 

9’2 

15-9 

E1 

X-1 

0 

1 

4 

9 

16 
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Solution  (3)  "Generalised  elliptical"  aaas  taper 


Blade  naas  per  unit  length 
Centrifugal  force  P  =  ^x  (l  -  x*)ndx 

=  2(n  7 1) 


»  -  (l-*2)n 

L  /,  _2\na 


1 -  (1  ..  x2)n  +  1 


The  generrl  differential  equation  (Hef.) 

&*£>  ■  -*£ 


becomes 

2GT7T 

<L| 

dx  | 

!<i-*2)“*1 

dxj 

O 

CM 

K 

l 

rH 

+ 

or 

1  -  X2  ( 

*2TToTS 

& 

=  0 

1'friting  5  for  (l  -  x)  [i.e.  transposing  origin  to  blade  tip] 

§£rf}6  * » -  »>  S  *  *  ■  ° 

The  solution  In  ascending  powers  of  x  is 

y  =  a,,  +  axx  +  .  .  .  .  +  + . 


Where 


ap 


s^LL-  (££3&j^  2 - k2} 

»(n  +  p;  n  +  1 
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Solutions  (A)  and  (B)  ■  ■  (l  -  x2)n 


2 

Deflection  y  «  Aq  +  a^ac  +  agX 

+  ...... 

+  ‘p^  + . 

(B) 

(A)  (n  =  i) 

2 

Where  »  -k  aQ 

*1  *  -k2‘o 

*2  =  (1-k2>  ^ 

a2  =  K>  t1  ~  fc2^  “l 

‘3  ■  M(M-“2) 

1  a2 

*3  -  $  <!  -  *2)  *2 

\  =  TOT  (if^f  _k 

)a3 

a4  =  ft  (5  “  fc2)  a3 

a5  =  (^Kt0 

3% 

5  (C  -  k  )  a^ 

a6  -  sfSsy  (^"1( 

0  ‘5 

a6  ■  i  <*?  ■  k2)  a5 

_  n+ 1  f§ru21  .  J. 
^  "  7Tn+tf  V  nJ.  51 

3** 

•?  =  35  (A6  -  k  )  ag 

„  n  +  1  /7h+28  , 

8  sfn+^J  ^  nd  k 

3*? 

a8  =  ife  -  k2)  *7 

.  _  .ALL  /'i£±2i  .  v 
a9  _  9Gw5T  V  n+1 

3* 

9  =  57  3  "  k  ^  a8 

n+1  /  9n+ii.5 

*10  x  lO(n-flO)  ^  ml 

*10  "  (33  -  k2)  a9 

i 
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EUlptioal  Taper  (Solution  A) 


gESSMSBC 


k2  si 


k2  -fi 

K  "  3 
k2=5 


k2  .  8 


k2  -2 
k  -  3 

k2  =  16 


Mode 


y  »  1  - 


y 

y 

y 

y 

y 


1  -  5x  ♦  fix2  -  2X3 
1  -  fix  +  16-ebc2  -  12-ax3  +  Vlx* 


1  -  f«  .  Jfi,2  -  MS3  .  -  ^c3 

1  .  16,  i  72«2  -  Sffc3  .  fW*  -  *  fx6 


The  above  are  natural  modes  (symmetrical  or  anti-symmetrical ) 
far  a  complete  two-bladed  rotor.  The  anti-aymmetrioal  modes 
(k2  =  1,  5,  ^  ,  .)  apply  also  to  a  single  blade  hinged  at  the 

»-h«  of  rotation  (x  si). 

For  other  hinge  positions  In  general  we  have  Intermediate 

o 

values  of  k  ,  leading  to  a  non-terminating  series: 
e.g.  k2  s  1*5 

y  »  1  -  l*5x  +  0-225x2  +  0-0375*3  +  0*01090bc^  +  0*00383^ 


♦  0-00152X6  +  0*00063x7 


This  is  found  to  have  a  hinge  (y  ■  0)  at  x  =  0*770 

i.e.  0.23R  from  the  axis 


I 


Pfcgs  U 


Elliptio al  Taper,  oase  (A) 


Vibration  amplitude  j  ail  function  of 

o 

frequanay  (k  )  and  bl*3s  station  5 


5 

6 

7 

12 

3 

14 

16 

1 

0 

+•123 

+•185 

0 

—118 

-143 

0-9 

-.098 

+•030 

+.124 

+•095 

-034 

—110 

0*0 

-184 

-•072 

+•029 

+  159 

+•066 

-017 

0-7 

-246 

-•167 

-.076 

+■169 

+•145 

+•080 

0*6 

—272 

—238 

-174 

+•113 

+•162 

+•150 

0*5 

1 

to 

VJl 

0 

—263 

—  242 

0 

+•097 

+•143 

0*4 

-168 

—225 

-249 

-137 

-039 

+•037 

0-3 

-014 

—103 

-166 

*0 

CM 

1 

-188 

1 

• 

F* 

0*2 

+•224 

+•126 

+•044 

-183 

-•22 4 

-233 

0*1 

+•558 

+•486 

+•413 

+•  l6l 

+•064 

-005 

0 

•a 

■a 

+1 

+1 

+1 

+1 

1 
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Solution  (D)  "Generalised  perebolto"  naas  taper 

Blade  nn  per  unit  length  n  =  (l  -  x)n 
Centrifugal  f owe  P  *  '\{l  -  x)n  dx 

'ft.-.)"-1  a  -jd1"8 

n  +  1  n  +  2 

The  general  differential  equation  (Ref. ) 

4  (pg)  - 

»—  4 (WS1  S ] * kV .  o 

Mw  p  2 

1  +  nr  -  n+1  x  4v  dY  ,_2 

®  TnfffCT -  3S?“aS  +  ^  =  0 

Writing  x  for  (l  -  x)  |i.e,  transposing  origin  to  blade  tip  ] 

x(sTT-rh)^  +  (1“x)S  +  k2y  *  0 

The  solution  in  ascending  powers  of  3  is 

y  =  »0  +  V  + . +  “p  *p  . . 

*p  "  (*'V".Va  *  "0  Vi 
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Parabolic  taper,  out  (C) 


Blade  vibration  aaplitude  2  as  a 
function  of  frequency  (k2)  and  blade  station  (x) 


.  k2 . 

4 

5.4 

7 

8*8 

11 

13 

15 

x"^- 

1 

-U9 

+•080 

+•182 

+•143 

+.004 

-.099 

-132 

0.9 

-198 

-006 

+•131 

+  •174 

+.096 

-Oil 

-092 

0-8 

-•250 

—114 

+•047 

+•149 

+•155 

+.085 

-008 

0-7 

-269 

-195 

-070 

+•077 

+  154 

+.348 

+•095 

0-6 

—248 

-246 

-157 

-•026 

+•095 

+.146 

+.151 

0-5 

-182 

-249 

-229 

-145 

-.026 

+.0 66 

♦•125 

0-4 

+•066 

-189 

—2 41 

—227 

-155 

-.071 

+.008 

0-3 

+•106 

-050 

—162 

—221 

-233 

—203 

-154 

0*2 

+•339 

+  184 

+•045 

—069 

-l6l 

-203 

-227 

0*1 

+•635 

+•528 

+•419 

+•308 

+•193 

+.104 

+•029 

0 

+1 

+1 

+1 

+L 

+1 

+1 

+1 
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Solution  (D)  extended 


The  mass  distribution  m  =  (l  -  x)n  Is  unsymmetrlcal 
about  the  origin  (x  =  0)  so  that  a  different  blade  taper  Is 
represented  by  the  expression  with  2  negative.  Convenient 
solutions  for  the  vibration  modes  oan  be  obtained  for  these  'blade 
extensions 1  also,  and  these  are  of  practical  value,  particularly 
for  n  <  0,  since  the  corresponding  taper  curve  is  hollow  and  has 
a  finite  tip  ordinate. 


The  length  of  the  extension  has  to  be  ohosen  to  sake 
F  =  0  at  the  new  blade  tip,  i.e.  the  origin  (x  =  0)  has  to  be  at 
the  o.g.  of  the  total  mass.  From  the  expression  for F  already 
derived,  the  length  of  the  extension  Is  seen  to  be  •— -r  , 

1  n  +  x 

i.e.  x  =  -  -j-  at  the  new  blade  tip. 

A  further  tramsposition  of  the  origin  to  this  point  is 
found  desirable  so  as  to  give  convergence  of  the  series  solution. 
The  new  differential  equation  is 

and  the  solution  in  ascending  powers  of  x  is 

y 


a0  +  V  + 


V?  + . 


a  = 
P 


p* 


V  n  +  2  *  J  Vl 


with 


Case  (D) ,  extended  (^3) 


The  "modified  parabolic"  mass  distribution  (l  -  x)n 
with  n  ranging  from  0  to  1  o overs  the  uniform  (n  ■  0)  end 
uniformly-tapered  blade  (n  =  l)  and  as  intermediate  oases  parabolas 
of  various  orders,  including  the  ocmmon  parabola  (n  = 


When  n  is  negative  the  mass  becomes  infinite  at  the 
blade  tip  (x  =  l)  but  the  solution  remains  valid  and  useful  for 
negative  values  of  x,  i. e.  for  the  blade  continuation  beyond  the 
cods.  The  special  value  of  this  result  is  that  the  mass 
distribution  curve  is  now  oonoave,  thus  more  resembling  aotusl 
blades,  and  further  it  has  a  finite  ordinate  at  the  blade  extension 
tip. 


The  latter  is  given  by  (l  -  x) 

x 


This  is  the  distance  from  the  axis  at  ■which  the 
centrifugal  force  vanishes; ;  in  other  words  the  axis  of  rotation 
is  at  the  centroid  of  the  blade  plus  its  extension.  The  tip 
ordinate  is 

The  taper  and  concavity  range  is  seen  from  the  table. 


Hass  curve 

t* 

II 

o 

-& 

-1 

(-D 

Boot  ordinate 

1 

1 

l 

(  ") 

Hid  ordinate 

1 

•707 

•438 

(2«0 

Tip  ordinate 

1 

•577 

•299 

<*) 

Area  coefficient 

1 

•732 

•495 

0 

Centroid  from  axis 

0-5 

•455 

•404 

0 

The  limiting  oase  (n  =  -l) ,  which  requires  to  be  solved 
separately,  is  also  quoted  for  comparison. 
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lyperbolio  nut  distribution  (E^) 


Blade  mass  per  unit  length  =  x-1  (x  from  aids) 

or  (l  -  x)"1  (x  frcei  tip) 

Tliis  oan  be  regarded  as  a  special  case  of  the  extended 
solution  (D) ,  with  n  =  -1,  but  this  gives  an  infinite  blade  length 
A  solution  far  finite  blade  length  (=  l)  is  however  readily 
obtained,  namely,  in  ascending  powers  of  distance  x  from  the  tip. 

y  =  aQ  +  a1x  +  ......  +  a^  + . 

*“•  ftp  =  jvi 

The  faot  that  the  mass  grading  becomes  infinite  at  the 
axis  means  that  the  length  of  blade  to  which  the  solution  oan  be 
applied  has  to  be  cut  short  of  this  point. 

Regarded  as  a  special  case  of  the  extended  solution  (D) 
it  is  seen  that  the  hyperbolic  mass  distribution  gives  the  limiting 
degree  of  hollowness  of  the  mass  curve  to  which  the  solution 
applies. 
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She  characteristic  bend  stiffness  end  shear  stiffness 


She  various  exaot  solutions  worked  out  for  rotating  flexible 
blades  apply  it  Is  found  also  to  non-rotating  blades  having  the  sane 
mass  distribution  end  a  certain  oharacterlstio  bend  stiffness 
distribution  related  to  the  mass  distribution.  Similarly  for  shear 
stiffness. 

It  is  considered  sufficient  to  uark  out  one  typioal  example. 


For  case  (D),  mass  distribution  (l  -  x)n  (x  from  root) 

and  for  the  first  non-rigid  mode,  with  A 

(applicable  to  the  blade  proper  plus  its  extension  beyond  the 
axis,  with  both  ends  free). 


The  node,  an  substitution  in  the  general  series  is 

y  *  i  -  x  +  & a  x2  (x  from  tip) 

The  shear  force  is  proportional  'o 

lyz?  a*  -  — i—  xn+1  _  ZiR-tll  (n  +  lj.fa  ±Jtl  x«+3 

“  -  n  +  1  x  (n  +  2)2  X  +  (n  +  Z)*  X 

By  further  integration,  the  herding  moment  is 


n-*2  2 

'  (n  +  W 


rnTiJ^’+TT 1  •Trrwxn+5+7^ 


.T.  rt...  T 

(n7WX 


n+if 


The  curvature  being  constant,  this  is  also  the  required  bend 
stiffness 


El 


-  3t^+2\n-+ *2  -  n  V*1  x)  apart  tr am  a  constant  multiplier. 


Similarly,  the  shear  stiffness  got  by  dividing  the  shear  force 
by  the  slope  is 
a+1  / 


GA.  =  x 


(—  -  n  V  T  xj  i'ron  a  constant  multiplier. 


Analagous  results  are  found  for  all  the  types  of  mass 
distribution  investigated,  and  these  ere  found  to  be  general  and 
applicable  to  all  modes. 
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Hinge  positions  versus 


Uniform  A  Tapered 


k2 

X 

k2 

X 

/ 1 

0.000 

f  1 

0-000 

\  1-1 

0-063 

1-3 

0-152 

Oth  / 

Oth/ 

0-167 

)  1*5 

/ 

0-231 

(1<5 

0-250 

V2 

0-3 U2 

f6 

0-000  ■ 

-0-012,. 

1th  .  7 

0-112 

1th  ^  5 

I 
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Application  of  results  to  fixed-root  blades 


The  results  for  anti-symmetrical  vibration  nodes  (y  =  0 
at  x  s  l)  of  a  rotor  oonsisting  of  two  similar  blades  (eg,  Cases 
A,  B)  and  for  the  corresponding  nodes  in  other  oases,  oan  be  regarded 
as  applying  equally  to  the  anti-synmetrio  vibration  when  the  two 
blades  are  rigidly  connected  to  each  other  via  a  rigid  hub,  subject 
to  possible  restraining  effects  on  the  hub  from  its  mounting  on  the 
vertical  driving  shaft  (i.e.  the  hub  may  not  be  freely  pivoted,  os 
assumed,  at  the  axis  of  rotation). 

However,  we  are  interested  also  in  applying  the  results, 
as  far  as  feasible ,  to  the  symmetrical  vibrations  of  such  rigid -hub 
type  rotors,  vjhen  the  root  condition  for  the  vibration  amplitude  or 
blade  deflection  is  ^  =  0.  If  the  hub  itself  is  sufficiently 
rigid,  this  root  condition  may  be  applied  at  the  station  -where  the 
blade  root  enters  the  hub. 

Such  a  condition,  in  conjunction  v/ith  y  =  0  (which  will 
be  true  if  the  blade  mass  is  small  enough  in  relation  to  that  of  the 
rest  of  the  helicopter,  and  the  rotor  mounting  rigid.)  oannot  be  met 
by  any  curve  applicable  to  a  perfectly  flexible  blade,  as  here 
assumed,  but  an  equivalent  hinge  position,  './hero  -the  condition  y  =  0 
cnly  has  to  be  met,  oan  be  determined,  making  use  of  the  foot  that 
the  restraining  effect  of  the  hub  on  the  blade  bending  extends  only 
over  a  limited  distance ,  which  oan  be  estimated.  It  is  in  fact 
Ve57?  very  nearly,  as  for  a  blade  of  uniform  stiffness  (El)  and 
constant  tension  (P);  this  is  the  distance  of  the  equivalent  hinge 
(or  effective  node)  from  the  point  where  the  blade  enters  the  hub. 
This  is  the  actual  distanoe,  to  be  expressed  as  a  fraction  of  the 
rotor  radius  R  in  accordance  with  the  terminology  here  used. 

Rigid  mounting  of  the  hub  in  the  vertical  direction  (y  =  0) 
has  been  assumed  but  if  this  condition  is  not  fulfilled  the  mode 
position  will  be  o hanged  accordingly,  by  an  amount  dependent  on  the 
sparing  oonstant  of  the  mounting  and  the  mass  of  the  hub  in  relation 
to  that  of  the  blade  In  general  the  effeotive  node  will  tend  to 
move  radially  outward  to  some  extent.  This  applies  for  instance 
-,/hen  the  mounting  is  flexible  enough  to  leave  the  hub  substantially 
free 
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